Helium in superstrong magnetic fields 



O.-A. Al-Hujaj and P. Schmelcher 
Theoretische Chemie, Institut fiir Physikalische Chemie der Universitdt Heidelberg, INF 229, 69120 Heidelberg, Germany 

(December 21, 2002) 

We investigate the helium atom embedded in a superstrong magnetic field 7 = 100-10000 an. All 
. . . effects due to the finite nuclear mass for vanishing pseudomomentum are taken into account. The 

' influence and the magnitude of the different finite mass effects are analyzed and discussed. Within 

»' \ our full configuration interaction approach calculations are performed for the magnetic quantum 

■ numbers M=0, — 1,— 2,— 3, singlet and triplet states, as well as positive and negative z parities. Up 

to six excited states for each symmetry are studied. With increasing field strength the number of 
, bound states decreases rapidly and we remain with a comparatively small number of bound states 

for 7 — 10'* an within the symmetries investigated here. 
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I. INTRODUCTION 



i-C ' The term "strong" field characterizes a situation for which the system is in the nonperturbative regime, i.e. where 
the magnetic forces are of the same order of magnitude or greater than the Coulomb binding force. For the ground 

a state of the hydrogen atom this corresponds to field strengths 7 > 1 a.u. (1 a.u. corresponds to 2.35 x lO'"^ T). We 
_ , refer to the term "superstrong" to indicate a field strength of 100 atomic units and more. 

^ • The motivation to study atoms and molecules in strong magnetic fields originates from several sources. Certainly 
^ \ the properties of atoms and molecules in strong magnetic fields are interesting from a pure theoretical point of view. 
c/2 ■ Due to the competition of the spherically symmetric Coulomb potential and the cylindrically symmetric magnetic field 
. ^ [ interaction we encounter a nonseparable, nonintegrable problem already for a one-electron system, i.e. the hydrogen 
. atom. Therefore it is necessary to develop new techniques to solve the Schrodinger equation in strong magnetic fields. 
The discovery of strong magnetic fields on the surface of magnetic white dwarfs (10^-10"* T) and neutron stars (lO**- 
10^ T) is a further major motivation. The spectra of these astrophysical objects can be dominated by the influence 
of magnetic and electric fields. For the analysis of atmospheres of magnetic white dwarfs and neutron stars, it is very 
important to have reliable data on the behavior of matter in strong magnetic fields. As an example we mention the 
the white dwarf GrW-f 70°8247. The interpretation of its spectrum was very important for the understanding of the 
properties of spectra of magnetic white dwarfs in general (see Refs. 

Highly accurate data are available for hydrogen in strong magnetic fields (see e.g. 1^,0 )■ This system is now 
understood to a very high degree. But there is also a significant interest in accurate data on heavier elements such 
1: as He, Na, Fe and molecules. Especially helium plays an important role in the atmosphere of certain magnetic white 
O ' dwarfs (see e.g. 

CN ; There were several attempts to calculate accurate energies for bound states of helium, but for astrophysical appli- 
• cations an accuracy of approximately lO^'' for the energies is needed for a large number of levels. We will concentrate 



c/3 \ here on investigations that address the high field regime: In 1975 Mueller et al. |11| calculated the few lowest levels 
i of He for 7 up to 20000 a.u. using a variational approach. Virtamo |l^ presented Hartree-Fock calculations on 
' the ground state (which is a triplet state with magnetic quantum number equal to —1). The same state has been 
considered by Proschl jl^ in 1982 in the range 21 - 21000 au. Vincke and Baye jl4| provide correlated calculations 
■ (7=4.2, 42 and 420 a.u.) for the lowest singlet and triplet states with positive z parity and magnetic quantum numbers 
rT'l M=0,— 1,— 2. In the work of Thurner [|l5| several triplet states are considered in the very broad range 7=8 x 10~^ — 
K*" ' 8 X 10^ a.u. We mention also the important work by Becken and Schmelcher, which covers the same symmetries 
k>( ] and the same number of excited states in the broad range 7 = 8 x 10"'' — 100 a.u. of the magnetic field strength 
• published in a series of papers [pj[-p^ . 

' The properties of matter in superstrong magnetic fields are especially interesting for the physics of cold neutron 
stars [p0|-p2t. But for the case of superstrong magnetic fields the finite nuclear mass effects become increasingly 
important. This is due to the fact, that energy shifts, caused by the finite nuclear mass, are of the order of 7/M0, 
where Mq is the mass of the nucleus. Evidently this correction is for superstrong magnetic fields of the same order of 
magnitude as the ionization energies. 

Of course the conceptual and in particular the computational situation becomes more complex when the full 
Hamiltonian, i.e. for finite nuclear mass, is taken into account. For both neutral (2^j2j] as well as charged systems 
p3|-p2| we encounter couplings among the different electrons as well as couplings between the collective and electronic 
motion. 
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Some approximations to the ionic problem are available: An approximate separation of the collective and relative 
motion of the charged system has been introduced in |2^ ] and applied to calculate the finite nuclear mass corrections 
for low- lying levels of hydrogenic ions in a magnetic field . Elaborated MCHF computations as well as adiabatic 
approximations were given in Ref. However these approaches cannot be used to calculate accurate results for 

all field strengths and all states investigated here. We mention that the behavior of the ion He+ has to be known in 
order to determine the ionization energies of the neutral He atom. 

In the present paper we provide a full configuration interaction calculation for helium in superstrong magnetic 
fields. All finite mass effects at zero pseudomomentum are taken into account and are analyzed. In section ^ we 
will describe the Hamiltonian, and some technical details co ncer ning our calculation. Furthermore we provide some 



remarks on the problem of the threshold energies. In section III we analyze the deviations of the Hamiltonian in the 



infinite nuclear mass frame from the full Hamiltonian. Ionization energies and transition wavelengths, resulting from 



our calculations, are provided in Sect. IV 



II. FORMULATION OF THE PROBLEM 

A. Hamiltonian and Symmetries 

To investigate the He atom we use a nonrelativistic approach. This is well justified by the fact, that relativistic 
corrections have been shown to be very small in strong and even superstrong magnetic fields For the sake 

of simplicity we take the electronic spin g factor to be 2, but our results can be easily adapted to any g factor 
by multiplying the spin operators and their eigenvalues by g/2. On the other hand, the ionization energies, and 
transition wavelengths, which are presented in this paper are not affected by this choice. The magnetic field vector 
will be denoted by B, whereas its magnitude will be denoted by 7 = \B\. The magnetic field vector B is chosen to 
point in z direction. 

The first step in ou r appr oach is the pseudoseparation of the collective and relative motion for the Hamiltonian in 
the laboratory frame ||2^|2^ , ^ which exploits the conservation of the so-called pseudomomentum K. The resulting 
transformed Hamiltonian is divided into three parts, which are denoted by and iJs. The operator Hi ~ 

K"^ /{2Ma) involves only center of mass (CM) degrees of freedom, where Ma is the mass of the atom. iJs contains 
exclusively electronic degrees of freedom. The operator H2 = K /Ma ■ B x represents the coupling between Hi 

and ffa, i.e. the CM and electronic degrees of freedom. It involves the motional electric field 1/Ma{B x K) which 
arises due to the motion of the (neutral) atom in the magnetic field and is oriented perpendicular to the magnetic 
field. This coupling is proportional to the pseudomomentum, and therefore vanishes for vanishing pseudomomentum. 
The pseudoseparation is possible for neutral systems only, since only then all components of the pseudomomentum 
commute. 

Within the present work we assume a vanishing pseudomomentum. Therefore we have no additional motional 
electric field and all effects due to the finite nuclear mass are included in H^ which is a function of the mass of the 
nucleus Mq and the magnetic field 7. In atomic units the electronic Hamiltonian takes the following form (internal 
coordinates are taken with respect to the nucleus): 



where 



H{Mo,-f)=H„n + H„,p (2.1) 



1=1 



1 



|ri - r2\ 

and 



(2.3) 



= (f-f. ■ X + . (2.4) 

The reduced masses are fi = 1/(1-1- 1/Mo) and fi' — 1/(1 — 1/Afo)- The Hamiltonian can be considered to consist of three 
parts. The first part contains the one-particle operators l/(2/i)p|, the Zeeman terms l/{2fi')B ■ li, the diamagnetic 
terms l/(8/i)(B x r^)^, the attractive Coulomb interaction with the nucleus — 2/|ri|, and the spin Zeeman terms 
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B ■ Si. The second part contains the two particle operator (2.3), which describes the repulsive Coulomb interaction 
between the two electrons. The third operator is the so-called mass polarization operator Hmp- It arises due to 
the transformation of the laboratory coor dina tes to the internal coordinates, which are relative to the nucleus. The 



reader should note that the Hamiltonian (2.1) has the same good quantum numbers as the Hamiltonian for infinite 
nuclear mass: the total spin 5'^, the z component Sz of the total spin, the magnetic quantum number M and the total 
spatial z parity II^ (parity is not an independent symmetry, it can be deduced from these corresponding symmetry 
operations). In the following we will denote the states by u^^^^M^" , where 25+1 is the spin multiplicity and 
= 1,2,3... denotes the degree of excitation within a given subspace. 



If we compare H{Mq,^) in Eq. (2.1) to the electronic Hamiltonian of helium in the infinite nuclear mass frame 
p6t i?(oo,7), we observe two different kinds of corrections due to the finite nuclear mass. The first is due to the 
occurrence of the reduced masses /i and n' in i?!-™ which are the so called normal mass corrections. The spectrum 
of the Hamiltonian Hrm which contains exclusively these normal mass corrections can be related to the spectrum of 
the Hamiltonian H{oo,"/') at a different field strength 7', via a unitarian transformation and an additional trivial 



energy shift (see section |I| and also Refs. [|6|,|16|,^ ) . The second type of finite nuclear mass effects is due to the 
mass polarization operators Hmp- We will call these specific finite mass corrections. They are by no means trivial 
and are related to the correlation of the electrons (the corresponding operators contain one-particle operators of both 
electrons). If the electrons behave in a correlated way, the specific nuclear mass effects are enhanced. Note that these 
operators are a consequence of the transformation from the laboratory to our non-incrtial system thereby eliminating 
the CM degrees of freedom. 

B. Technical remarks 

Some comments concerning our computational approach are in order. Its basic ingredient is an anisotropic Gaussian 
basis set, which was put forward by Schmelcher and Cederbaum pTj . This one-particle basis is sufhciently flexible to 
to describe finite electronic systems for any field strength and was successfully applied to several atoms and molecules 
in magnetic fields [^-18 3S 



This basis set is not only flexible: in the case of atoms all matrix element can be calculated analytically and in par- 
ticular evaluated efficiently ||l^,|l^. The price, which has to be paid, is that for each field strength and each symmetry, 
the basis set has to be (nonlinearly) optimized. This is achieved by variationally computing the eigenfunctions of the 
corresponding one-particle problems (H, He+, etc). This task is time consuming and needs experience, because the 
starting values for the nonlinear variational parameters have to be chosen carefully in oder to obtain well-converged 
results. In recent times however this has become more comfortable and was speeded up by new optimization algo- 
rithms and tools, which allow an almost automatic construction of a basis set. The implementation of these tools 
still needs a proper selection of the coefficients and an evaluation of the outcome of the optimization. The latter is 
of particular importance for the field regime addressed in the present work for which small changes in the basis set 
can hav e a considerable effect on the results of the full problem, i.e. the eigenenergies and eigenfunctions of H{Mq, 7) 



in Eq. (2.1). A one particle basis set typically consists of 200-400 basis functions, which combine 3000 to 5000 two 
particle configurations. The full CI approach leads then to a generalized eigenvalue problem. Numerical problems 
arise in this generalized eigenvalue problem if near linear dependencies in the basis set arise. These dependencies 
cannot fully be avoided, but the resulting numerical instabilities can be removed by a cutoff of the small eigenvalues 
of the corresponding overlap matrix. 

The most CPU time consuming part in the above approach, is the evaluation of the electron-electron matrix 
elements. Although analytical formulas for all matrix elements are available, it is important to have efficient algorithms 
for their evaluation, since the matrix elements for the electron-electron Coulomb interaction are by no means trivial. 
A detailed and sophisticated analysis of their analytical representation is crucial. In the simplest form it contains 
multiple sums over hypergeometric functions Jl6| , pTt . For the evaluation of the hypergeometric function all possible 
analytical continuation formulas have been worked out (see Ref. p^) and made it possible to reduce the time for the 
calculation of the matrix elements by a factor of 50 compared to a straightforward implementation. We emphasize 
that the present investigation would have been impossible without this efficient implementation of the evaluation of 
the matrix elements. 

C. Threshold 

The result of the solution of the generalized eigenvalue problem are the total energies of the system. These energies 
are however not of primary interest: they increase almost linearly with the magnetic field strength due to the raise 
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of the kinetic energy in the presence of the field. The relevant energies are e.g. the transition and ionization energies 
which are for the high field regime addressed here by several orders of magnitude smaller than the total energies and 
therefore 'hidden' in the total energies. To obtain the ionization energies we need to know the energy of the He"*" 
ion taking into account the finite nuclear mass in the magnetic field. These energies for He"*" are unfortunately not 
known exactly. In contrast to this the energies for He"*" assuming an infinite nuclear mass can be calculated from 
the well-known energies for the hydrogen atom with fixed nucleus in strong magnetic fields (see e.g. Ref. 0) via the 
corresponding scaling relations (see for example Q). However the energies of the He"*" ion with finite nuclear mass 
cannot be extracted from results for the hydrogen atom (with finite or infinite nuclear mass) because of the unique 
coupling between the collective and electronic motion which has to be taken into account. This coupling inherently 
mixes the collective and electronic motion and the exact ground state therefore combines both motions. The moving 
He"*" -ion in a strong magnetic field is a problem which has not been solved numerically exactly in the literature. As 
a result the exact threshold energy for He"'' + e~ is not known. Using the threshold energy of He"*" + c^ for fixed 
nucleus E^^'j^ provides a wrong description of the threshold energies in superstrong fields. 
The Hamiltonian for the helium positive ion reads in atomic units as follows: 

Ho = Ha+Hb + (2.5) 
Ha = TTTIT^ (P-^B^ r) (2.6) 



2{Mo + 1) V 2 

12,1/, \ „ , 1 / ,2 , 



Here Ha describes the collective motion of the ion as a free particle with charge 1 and mass Mq + I moving in a magnetic 
field {R and P are the center of mass coordinate and momentum, respectively). The operator Hb couples electronic 
and collective motion. The operator He is the electronic part of the He+ Hamiltonian with cji = 1 + 1/{Mq + 1)^ and 
602 = 1 + (2Mo + l)/(A/o + 1)2. 

Although the exact threshold energy for He+ + e~ is not available some approximations to it have been calculated 
in the literature. One of those approximations E'^^^ ignores the coupling between the collective and electronic motion 
Hi,. The corresponding threshold values consist of the sum of eigenenergies of the electronic Hamiltonian He and the 
zero point energy for Ha- The zero point energy for the collective motion is the corresponding energy for the lowest 
Landau level. The energy for this threshold is typically too high, i.e. the true ionization energies are overestimated, 
since the coupling Hh, which is neglected in this approximation, tends to reduce the threshold energy. All ionization 
energies shown in the present work are calculated by applying E*p^^. 

The second approximation we will choose, ignores the zero point energy of the ion and the coupling term, i.e. Ha 
and Hb- Therefore, only the eigenenergies of the mass corrected electronic Hamiltonian He are taken into account. 
This threshold is denoted as E^^. It is motivated by the fact, that for an infinitely strong magnetic field the energetic 
contributions due to the zero point energy of Ha and the coupling Hb exactly cancel. A third alternative threshold 
energy E^^ is obtained by employing the adiabatic expansion approach presented in |3^ . Since this approximation 
only takes into account the lowest Landau energy for the CM motion, it becomes increasingly accurate for increasing 
energetic separation between the Landau levels (and therefore increasing magnetic field strength). For field strengths 
below 2000 a.u. this approximation is not reliable. It can be observed, that the threshold E^^ approaches the threshold 
energy E^^^ in the limit of high fields. E^'^ is expected to be very accurate for sufficient high field strengths. 

III. THE FINITE NUCLEAR MASS EFFECTS 

There is a transformation, which connects the spectrum of the infinite nuclear mass Hamiltonian H{oo, 7) with the 
spectrum of the Hamiltonian H^m fl^ ]. This transformation reads as follows: 

Here U denotes a unit arian transformation, which transforms r ^ r j [i and p — > p/i. The second term on the right 



hand side of Eq. ( |3.l| ) represents a field dependent trivial energy shift, since the total z component of the spin and 



the total z component of the orbital angular momentum are conserved quantities. 
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From equation ( |3.lD it can be seen, that the leading mass correction is of the order of 7/M0. For the states with 
magnetic quantum number M < 0, this correction is positive. This means that the corresponding energy levels 
are shifted i.e. increase linearly with 7 and eventually pass the ionization threshold. The latter implies that the 
corresponding bound electronic state becomes ionized. This shift depends exclusively on the magnetic quantum 
number M, the total z component of the spin S^, the magnetic field strength 7, and the nuclear mass Mq. 

The energy difference A_Erm=Erm(7)"E(oo, 7), where Erm denotes the eigenvalues of Hrm and E(^oo, 7) the eigen- 
values of i?(oo,7), can be related to the energy difference Ai?^„j of the eigenvalues of the Hamiltonian H^- and HJ;,„. 
The operator Hj describes two free noninteracting electrons 



Hf-T. \p\ + + ^P? (3.2) 

i—l 

whereas refers to the corresponding 'artificial' Hamiltonian with reduced masses 



2 2 
^^--E^^'^ + ^^^ + fc^?- (3-3) 
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For the energetically lowest Landau level with negative magnetic quantum number M and vanishing momentum in 
z direction we have 

^Krn = liT (1 + 1^1) ■ (3.4) 

It cannot be expected, that the normal mass corrections of the helium atom Ai?J;„j follow exactly equation {\i.4), 
because the Hamiltonian Urm contains the interaction between the two electrons and of the electrons with the 
nucleus. However it is suggestive to introduce a parameter 5(7) such that 

AErm = ErM " ^(oo, 7) = ^E';.„^{^) {1 + <5(7)} . 

For the states and field strengths investigated in the present work we have (5(7) ^ 1. Therefore 5 can be treated 
as a correction term. Since this correction is due to the Coulomb interaction it is state dependent. In Fig. |l| the 
quantity 5 is shown as a function of the magnetic field strength 7 for a few selected singlet states belonging to several 
different symmetries. It can be seen, that <5(7), for all states considered, follows a power law (^(7) = C ■ 7^'^, where 
the exponent A « 0.62 does not depend on the state. However the corresponding proportionality constant C varies 

over nearly one order of magnitude for the different states. 

To understand more on the behavior of the quantity 5 we expand the first part on the right hand side of Eq. (3.1) 
in powers of I/Mq. Omitting the spin part we obtain: 

Mi^(oo, ^) ^ i.™(7) + f ^ . {i^(oo, 7) + ^^^^ + 0(^)} . (3.5) 

_ g(oo,7) , ^ ig'(oo,7) 

-i^(oo,7)-^^ + 27^^. (3.6) 

Here the prime indicates the derivative with respect to the field strength. Now AErm and consequently S can be 
expressed in terms of the eigenenergies of H(oo, 7): 

-^(oo,7) + 27^^(oo,7)+7|M| 

WTm ' ^'-'^ 

Now the state dependency of d can be understood as the dependence on the derivative of the eigenenergies E(o o, 7) 



with respect to the field strength. We emphasize, that the quantities £'(00,7) and jE'{oo,j) in Eq. (3/7) are 
almost equal and therefore approximately cancel in the superstrong field regime. Thus ^(7) can be approximated by 
(5(7) « (£;'(oo,7) + |M|)/(1 + |M|) - 1. 

Fig. H (a) illustrates the mass polarization energies \Emp\ — \E{Mo, 7) — Erm{"t) \ for the energetically lowest singlet 
states and Fig. ^ (b) for the corresponding triplet states. Here E{Mo, 7) denotes the eigenenergies of H{Mq, 7). First 
we observe, that the absolute values of E^p are very small: For 7 = 10* they are at least eight orders of magnitude 
smaller than the corresponding total energies. They are also small compared to the normal finite mass corrections 
AErm- For 7 = 10''' Emp is at least four orders of magnitude smaller the normal finite mass corrections. Opposite to the 
quantity AErm the behavior of \Emp\ depends strongly on the state. For the states 1-^0"'",1^(— 1)+,1^(— 2)+,l^(— 3)+ 
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the quantity \Emp\ increases. These states contain so-called magnetically tightly bound orbitals ^,^,^| and therefore 
their ionization energy diverges logarithmically in the limit 7 — > cx). An increase of the influence of Emp\ can also 
be observed in Fig. |2| (b) for the states + 2)+ and 3) + , which represent also magnetically tightly 

bound states. For the remaining states 1^0^, 1)^ and the corresponding triplet states as well for 1*^0+ \Emp\ 
remains almost constant as a function of the magnetic field strength. This effect can be easily understood: For the 
magnetically tightly bound states the electrons are close to each other in a relatively narrow region of space and 
therefore electron correlation is important. The mass polarization operator is sensitive to electronic correlation, due 
to the fact that it contains operators of both electrons. The sign of E^p is not shown in Fig. ||. It is positive for 
the states related to the above mentioned tightly bound states and negative for the others. The only exception is the 
state l^O"*", which does not belong to the tightly bound states, but nevertheless E^p has a positive sign. 



IV. RESULTS 



In the following we will present our results for the ionization energies and transition wavelengths of the helium 
atom for magnetic fields ranging from 100 au to 10000 au. These investigations have been performed for the magnetic 
quantum numbers M — 0,-1,— 2,— 3, singlet and triplet states and positive and negative z parity. Only for the 
magnetic quantum number M = —3 exclusively positive z parity states have been studied. For most symmetry 
subspaces we investigated 6 excited states. 



A. Ionization energies 



According to the above the reader should keep in mind, that the exact energy of the ionization threshold is not 
known, and therefore the ionization energies are not known accurately. However the ionization energies calculated by 
using different approximative threshold energies Ej'^, E^^, El^^^, E^'^ show the same trend: the number of bound 
states of the helium atom becomes finite for superstrong magnetic fields in contrast to the situation without a magnetic 
field, or in the limit of an infinitely heavy nucleus, where the helium atom has an infinite number of bound states. Only 
the so-called magnetically tightly bound states are bound within the complete regime of field strengths and quantum 
numbers considered in the present work. We consider in the following the quantity Eion = E^J^^^ — E{Mo,j) as a 
function of the field strength together with the above mentioned approximations for the threshold. 

Fig. H (a) shows Eion for the six energetically lowest states of zero magnetic quantum number and positive z parity. 
The 1^0^ state is the most tightly bound state. In strong magnetic fields it represents however not the ground state 
of the atom, because energetically low-lying states are fully spin-polarized in the high field regime. The ground state 
is given by the 1^(— 1)+ state, which is also a tightly bound state. Fig. ^ (a) shows, that E^on increases for the tightly 
bound state 1^0"*", but remains approximately constant for all other states as a function of the magnetic field strength. 
Furthermore we observe that all states v^O^ with 1/ > 1 as well as the corresponding triplet states pass the threshold 
energies being either E*J^ or E'^^ with increasing field strength. The only remaining bound state for 7 = 10^ au is 
the 1^0+ state. 

In Fig. H (b) the corresponding quantity for states with M = and negative z parity are shown. For these states 
Eion as a function of 7 varies only to a very minor extent, which is due to the fact, that none of these states is a 
tightly bound one, and none of these states remains bound when 7 approaches 10"* au. For the triplet states with 
negative z parity the quantity Eion decreases slightly, for 7 > 100 au. This is not due to the finite mass effect, but can 
be also observed for the quantity Ej'^ — E{oo,^), whereas this quantity increases monotonically for all other states 
investigated in the present work. The reason is the complicated interplay between correlation, which tends to increase 
ionization energies and Coulomb repulsion, which tends to decrease it. For the states with M = the electrons are 
confined in a very small domain of space, which increases correlation as well as the Coulomb repulsion. On the other 
hand for triplet states the electrons are separated, because the wavefunction is antisymmetric, which reduces both 
effects. For the i/^0~ states the increase of correlation energy is smaller than the increase of the Coulomb repulsion 
energy. 

In figure ^ (a) Eion is shown for the states 1)+. For the magnetically tightly bound states 1^(— 1)+ and 

l'^(— 1)+ the ionization energy remains positive, i.e. these states are bound in the complete regime 7 < 10^ au. For 
higher excited states, i.e. 1)+ with v > 1 the energy E{Mo,j) becomes even larg er than E^J^^^ and therefore 

Eion decreases strongly on the logarithmic scale. To understand this we review Eq. (pTl): The dominant term on the 
right hand side is of the form — M7/M0 (the spin part does not affect the ionization energies). Therefore E{Mq, 7) for 
states with M < raises about this amount and will pass the threshold at lower field strength than their counterparts 
with M = 0. 
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In Fig. H (b) we show our results for Eion for M = —1 and negative z parity. A similar behavior to that of the states 
j/2S+i^_-[^-j + ^ with > 1 is observed: The ionization energy as a function of the magnetic field strength decreases 
rapidly. This is due to the fact, as mentioned above, that the finite mass corrections force these states to pass even 
the threshold and the states become unbound. For negative z parity tightly bound states do not exist. 

In Fig. @ (a) we present our results for Eion for the energetically lowest singlet and triplet states with M — —2 and 
positive z parity. Similar to Fig. ^ (a) the quantity Eion rises for the two energetically lowest singlet and triplet states 
from 1.8 au to approximately 10 au. These two states belong to the magnetically tightly bound states. However in 
contrast to Eion for the 2^"^+^(— 1)+ states in Fig. ^ (a) Eion for the state 2^(— 2)+, which is the first excited singlet 
state of this symmetry, does also rise and stays bound within the complete regime of field strength considered here. 
This is a remarkable feature since the influence of the finite mass effects for this state is even bigger than for the 
states with M = —1. The reason for this behavior lies in the presence of an avoided crossing which takes place at 
7 « 100 au. It can be seen that Eion for the higher excited singlet states is raised as well and it approaches the 
corresponding value for the next energetically higher triplet state. 

The ionization energies for the negative z parity states for AI = —2 in Fig. || (b) look similar to these of the M = —1 
in Fig. H (b). Magnetically tightly bound states do not exist for M — —2 and negative parity, therefore all states of 
this symmetry become unbound with increasing field strength. As mentioned above the influence of the finite nuclear 
mass increases with increasing magnetic quantum number \M\, therefore the states 2)^ become unbound at 

lower field strengths than the corresponding states with magnetic quantum number M = —1. 

In Fig. II we encounter that Eion for the energetically lowest singlet and triplet states with M — —3 and positive 
z parity is positive for all field strengths, considered in the present work. These two states belong to the magnetically 
tightly bound states. Similar to Fig. || (a) an avoid crossing takes place at 7 « 200 au. This can be more clearly seen 
if the quantity Ej'^ — E{oo, 7) is considered. But unlike to the case of M = —2, where Eio„ for the first excited singlet 
state is raised, here the triplet state 2'^(— 3)+ is shifted to lower energies and remains therefore bound for relatively 
high field strengths. Nevertheless, the state 2'^(— 3)+ passes the ionization threshold for 7 approaching 10^ au and 
therefore becomes unbound due to the infiuence of the finite nuclear mass effects. 



B. Transition wavelengths 

In contrast to the ionization energies, which strongly depend on the exact values for the threshold energies, transition 
wavelengths can be calculated from total energies without the knowledge of the exact threshold energy. The only 
feature which remains unknown for the transition wavelengths is the exact field strength, foe which the particular 
bound-bound transition disappears. Therefore we refer the total energy of a state to the threshold E^Jl^^ in ord er to 



decide upon its bound character as a function of the field strength. According to the discussion provided in Sect. II C 



the true field strength for which the state becomes unbound is lower than the value obtained by refering the total 

^ zpmc ■ 



energies to 



• • • • P I — I 

Some general remarks on our results for the transition wavelengths presented in Figs. |7| - |1J are in order. Linearly 

polarized transitions show the general feature, that there are two separated parts of the spectrum. A few transition 
wavelengths become shorter, following approximately a power law as a function of the field strength. Other transition 
wavelengths, which are much longer, stay almost constant as a function of the magnetic field strength. The short 
wavelengths correspond to transitions involving the magnetically tightly bound states. The other lines correspond to 
transitions between higher excited states. 

For circular polarized transitions the appearance of the spectrum of wavelengths is different. Only for circular 
polarized transitions which involve states with positive z parity a separated bundle of short wavelengths is present. 
This is again due to the presence of tightly bound states for positive z parity but their absence for negative z parity 
states. Transitions between higher excited states which do not involve tightly bound states, look different than their 
corresponding linear polarized transitions: the wavelengths increase or decrease as a function of 7. In the case of 
circular polarized transitions, states belonging to different magnetic quantum numbers M are involved and therefore 
the energy of the state with the higher absolute value of the magnetic quantum number M raises more strongly with 
increasing field strength than the energy of the state with the lower magnetic quantum number. This again goes 



back to the nuclear mass effects contained in Eq. (3.1) and yields the increase or decrease of the wavelengths for the 
circular polarized transitions with increasing field strengths. 

In Fig. the transition wavelengths for the singlet and triplet transitions among the j/^'^+^0+ and /i^'^+^0^ states 
are shown. The separation mentioned above, which is a general feature for linearly polarized transitions is obvious. 
The short wavelengths are the transitions involving the most tightly bound state 1^0+. According to Fig. ^ (a) its 
ionization energy increases with increasing field strength, whereas the ionization energy of all other (excited) states 
varies only marginally. 



7 



The circular polarized transition wavelengths shown in Fig. [s] involve the u'^^'^^O'^ and symmetry 
subspaces. Separated by a large energetically gap there is a bundle of short wavelengths, which is due to transitions 
to the magnetically tightly bound states 1^0+,1^(— 1)+ and 1)+ and a long wavelength part. 

The spectrum of circular polarized transitions among the subspaces j/^'^+^0~ and ^^'^"'"^(—1)", shown in Fig. ^jshows 
also the characteristics of circular polarized transitions discussed above. Since in this case only states with negative 
z parity are involved, we have no tightly bound states, and the part of the spectrum with very short wavelengths is 
missing. 

The behavior of the wavelengths of the linear polarized transition from the 1)+ states to the fi^^^^{—l)^ 

states is shown in Fig. |l^ and looks similar to the one given in Fig. 0. Two well separated parts of the spectrum can 
be distinguished: Short wavelengths, which are due to transitions to the magnetically tightly bound states 1)+ 
and l)"*" and decrease as a function of the field strength. Transition wavelengths between higher excited states, 
with wavelengths larger than 1000 Angstr0m, remain approximately constant. Due to the fact, that the influence of 
the finite nuclear mass is more significant than for those shown in Fig. ^, much less lines belonging to bound-bound 
transitions are present. 

The spectrum shown in Fig. |ll| for the symmetry subspaces 1)+ and /i^"^+^(--2)+ differs in some respect 

from the common pattern of circularly polarized transitions. Similar to other spectra circular polarized transitions, 
wavelengths involving transitions to the tightly bound states 2) + ,l'^(— 2) + ) can be easily iden- 

tified. As a function of the magnetic field strength, these wavelengths follow approximatively a power law. On the 
other hand there are structures shown in Fig. |ll|, which arise due to the influence of the finite nuclear mass. However 
in the gap between these two parts of the spectrum, additional lines occur which are caused by avoided crossings. 
The transitions i/'^^'^^ to iJ?^~^^{—2)~ shown in Fig. |lj show the clear signature described above for circular 
polarized transitions. Only a few lines belong to bound-bound transitions. 

The spectral transitions given in Fig. |l3| 2)+ to iJ?^~^^{—2)~) show the typical behavior of linear polarized 

transitions. Deviating from the general pattern, transition wavelengths to the state 2^ (—2)+ form their own bundle 
of short wavelengths, being located between the wavelengths of the transitions of the tightly bound states 2) + 
and 2)+ and the transitions among higher excited states. 

The spectrum of transitions among the subspaces 2)+ to 3)+ is shown in Fig. 
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Transitions 

to the magnetically tightly bound states 2)+,l'^(— 2)+, 3)+, 3)+ can be clearly identified. Also the 
transitions dominated by the normal finite mass effects can be seen for long wavelengths. Additionally the influence 
of the avoided crossings is visible. 



V. BRIEF SUMMARY 



We have presented the first systematic full CI calculations for helium in superstrong magnetic fields, taking into 
account the effects of finite nuclear mass. These effects are extremely important in the superstrong field regime, 
because the relevant parameter for the importance of the finite nuclear mass effects is j/Mq. We analyzed the 
influence of the normal and the specific finite nuclear mass effects. It has been shown that the leading finite nuclear 
mass effect, does not depend on the state, but only on the magnetic quantum number. The state dependent part of 
the normal finite mass effects depends on the derivative with respect to the magnetic field strength of the total energy 
of the corresponding state in the infinite nuclear mass frame. Furthermore it has been shown, that the specific mass 
effects, which are caused by the mass polarization operators, are very small compared to the total energies and small 
compared to the leading normal mass effects AE^^^. 

In the superstrong magnetic field regime, the spectrum of helium is cut off by the effects of the finite nuclear mass. 
We found that only a comparatively small number of 

states is bound in the complete regime of magnetic field strengths investigated in the present work. Although the 
exact ionization threshold for helium is unknown, all available approximations to the exact threshold confirm this 
trend. Transition wavelengths for many linear and circular polarized transitions were provided. Their typical behavior 
has been identified and the effects of the finite mass on the transition wavelengths has been analyzed. 

The determination of the critical field strengths (i.e. the field strengths where the individual states become unbound) 
require a detailed investigation of the ground state of the moving helium positive ion in a magnetic field. 
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FIG. 1. The parameter 5 as a function of the field strength for singlet states of various symmetries. It reflects the state 
dependent normal flnite mass corrections due to the Coulomb interaction. For details see text. 




FIG. 2. The absolute values of the mass polarization energies \E.mp\=\E{Mo,'y) — Erm{'y)\ for selected singlet (a) and triplet 
(b) states as a function of the magnetic field strength 7. 
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FIG. 3. Ionization energies Eion = -Bzpmc — E{Mo,j) of the helium atom for the jy^^^^Q^' states, (a) The ionization energies 
for the states with positive z parity, (b) ionization energies for the corresponding states with negative z parity. Note that the 
threshold E*J^^^ overestimates the ionization energies. Different approximations to the exact threshold E*^^, B*^ and the fixed 
nucleus ionization threshold are given. 
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FIG. 4. Ionization energies Eio„ = Elp^^ — E{Mq, 7) of the helium atom for the electronic states with M = — 1 as a function 
of the magnetic field strength 7. (a) Eion for positive z parity states, (b) Eion for negative z parity states. 
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FIG. 5. Eion for the electronic states v^^^^ {—2)^' . (a) For positive z parity wc encounter an avoided crossing, which gives 
rise to a significant increase of Eion for the 2^ (—2)"'" state, (b) Eion for the corresponding negative z parity states. 
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FIG. 6. Eion for the electronic states with magnetic quantum number M = — 3 and positive z parity. An avoided crossing 
increases the ionization energy for the state 2^ (—3)"*". 
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Fig. 7 
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FIG. 7. Wavelengths of linearly polarized transitions among the singlet and triplet states f^^'^^O'*' and IJ,^^~*'^0~ in Angstr0m 
as a function of the magnetic field strength in atomic units. The short wavelengths correspond to transitions involving the 
most tightly bound state l^O"*". 
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FIG. 8. Transition wavelengths for the singlet and triplet transitions from ^'^'^"'"^0"'" to 1)"*" in Angstr0m corresponding 

to circular polarized transitions. The shortest wavelengths are due to transitions involving the states 1^0'*',1^(— l)"*" and !''(— 1)^, 
which belong to the magnetically tightly bound states. The transitions between excited states are dominated by finite nuclear 
mass effects. 
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